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I. INTRODUCTION 


The investigation of thermoacoustic heat engines was begun by 
Wheatley and associates in the early 1980's [References 1 & 2]. 
They explored the phenomena upon which these engines are based by 
introducing a thermoacoustic couple (TAC) into an acoustic 
standing wave. A TAC is a stack of short, thin plates. At least 
one of those plates is instrumented with a series of 
thermocouples to measure temperature difference developed across 
it as a result of the heat transport. In 1987, Muzzerall [Ref. 
3] made measurements of the temperature difference developed 
across a TAC as a function of the position of the TAC in the 
acoustic field and compared the results to the theoretical value 
Given by Wheatley, et al. There was good agreement between the 
measured data and the theory at low acoustic pressure amplitude 
however, at high pressure amplitude there was agreement only in 
the regions near the pressure antinodes. Elsewhere, 
unaccountable severe skewing and an overall reduction in the 


temperature difference was observed. 


Muzzerall published extensive data only at two acoustic 
pressure amplitudes separated by 12 dB. Also, his spatial 
resolution was only 1 cm. The purpose of the experiment reported 
in this thesis is to extend this investigation to see where the 


theory and measured data start to diverge and to further explore 


what happens at high acoustic and mean pressures. In order to 
accomplish these goals, better spatial resolution and the ability 
to investigate a broader range of parameters is required. The 
parameters to be controlled are acoustic pressure amplitude, 
ambient (or mean) pressure, frequency, and gas type. In order to 
handle the quantity of data recorded and to ensure repeatability 
of measurements, both the data acquisition and the linear 
positioning of the TAC are computer controlled in this 
investigation. Following a derivation of the theoretical basis 
for the generation of temperature gradients and a summary of 
Muzzerall's results, the computer controlled experimental 
apparatus is delineated. Methodology for data acquisition is 
discussed and a graphical comparison of measured and theoretical 


data is made. Areas for further investigation are included. 


II. THEORY 


Wheatley, et al. [Ref. 2] used the term "natural engine" to 
describe a process where a phase shift is introduced between 
temperature and velocity by a natural, irreversible process - 
here thermal conduction. This process leads to, what has been 
termed, thermoacoustic heat transport. The TAC acts as such an 
engine establishing a thermal lag which drives an entropy flow. 
A brief description of this phenomena is described below. A 
detailed development of theoretical thermoacoustics applied to 


engines is given by Swift [Ref. 4]. 


Following the description of the phenomenon, a summary of 


Muzzerall's findings is presented. 


A. BASIC THEORY 


The cyclic heat pumping action of the TAC as a "natural 
engine" may best be described from a Lagrangian point of view. 
Referring to Figure 1, a poor thermally conducting plate is 
located near the end of a rigid resonant tube. A sound source is 
located at the left end of the tube (not shown). This source 
generates an acoustic standing wave in the tube. Consider a 
parcel of gas located near (within a thermal penetration depth) 


the plate. Both the plate and the parcel of gas are initially at 


Rigid Closed Tube 





Figure l 


The Thermoacoustic Effect 


temperature T. During the compression phase of the acoustic 
cycle the parcel is displaced toward the pressure antinode and 
undergoes an adiabatic compression increasing it's temperature to 
T++. The temperature difference drives a heat flow (Q) from the 
gas to the plate reducing the parcel's temperature to T+. On the 
expansion phase of the acoustic cycle, the gas undergoes an 
adiabatic expansion, but is now at the lower temperature T-. 

Here the temperature difference results in a heat flow in the 
opposite direction, returning the gas parcel to its original 
temperature T. The result is that after one cycle, the parcel 
has transported an amount of heat Q over a short distance. The 
deposited heat is picked up by an adjacent parcel and moved still 
further during the next cycle. The heat is shuttled down the 
plate, toward the closed end of the tube, until it reaches the 
end of the plate. There it collects, until conducted back 


through the plate or the surrounding gas. 


The system reaches steady state when the hydrodynamic heat 
(entropy) flow in the gas is balanced by the return diffusive 
heat (entropy) flow in the plate and in the surrounding gas. The 
entropy flow increases when there is no average temperature 
gradient in the plate. As the gradient increases, the heat flow 
begins to drop. There exists a critical gradient beyond which 
acoustics alone will no longer increase the temperature gradient 


in the plate [Ref. 1]. The gradients from zero through the 


critical value demark the heat pump regime, which is addressed 


specifically in this thesis. 


The behavior of the temperature difference as a function of 
its position in the standing wave can be derived from Wheatley's 
estimate for the time-average hydrodynamic heat flow [Ref. 1: p. 


Loo, Banzeao 1. 


o=r§, (7, 8)(7'-1)p,u, 
where 


7{ = overall perimeter of the plate 


b.. 


thermal penetration depth 


T, = mean temperature 

LB = isobaric expansion coefficient 
P, = acoustic pressure magnitude 

u, = acoustic velocity magnitude 


[* = defined by Wheatley as the ratio of Vr/V 0... 
The terms rB determine the adiabatic temperature change per 
pressure change and TO, is the volume of gas-plate interaction. 
Ce 1) shows how the entropy flow changes as the temperature 
gradient changes and becomes zero at Ver = Ve The next 


step in the development [Ref. 1: p. 158, Eqn. 30] is to note 


that: 


<Q, >/(KA/1) 


Ar 

where 
k = thermal conductivity 

A = sectional area of the plate 

1 = effective length of couple 


See) — Spatial average with subscript 2 denoting second order 
Wal Lael ty 


The relationship between acoustic pressure, velocity and the 


temperature difference is seen in Figure 2 for arbitrary units. 


Wheatley's expression for the temperature difference across 


fmemrAc [Ref. 2: p. 157, Eqn. 17} is: 
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OG = Prandtl # 
0. = thermal penetration depth 
a = sound speed 


K = wave number 
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Figure 2 
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Spatial Dependence of Pressure, Veloc 


X = distance from closed end of tube 

Po = mean density 

k ,kK. = thermal conductivity for plate(p) and gas(g) 

1 = length of TAC 

ae = TAC temperature 

G) = 271 . frequency 

<= ratio of specific heats, C,/C, 

which is what was used to determine the theoretical values in 


this paper. 
B. MUZZERALL'S RESULTS 


Included for comparison are Muzzerall's results (Figures 3 
and 4) for a TAC made of stainless steel plates 2.26 cm long x 
2.16 cm wide, at four penetration depth spacing. These 
measurements were made in argon at a frequency of 705Hz and a 
static pressure of one atmosphere. The lines indicated as theory 
are the predictions of the expression given previously for Ar. 
The acoustic pressure amplitude in Figure 3 is 630 kPa. There is 
overall good agreement between theory and experiment. In Figure 
4 the acoustic pressure amplitude has been increased to 2500 kPa. 
The overall agreement between theory and experiment is severely 
diminished, although it is still quite good in regions near 
acoustic velocity nodes. The impetus behind the present thesis 
is to explore the regions between and beyond these two figures in 


order to better understand the reasons for this behavior. 


~ MEASURED; —-- THEORY 
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Figure 3 
Muzzerall's Results 150dB 
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IIl. APPARATUS 


The experimental apparatus can most eaSily be described by 
breaking it down into its functional component parts. The 
subsystems are: the tube and pressure vessel, the TAC; the 
linear positioning system; and the data acquisition (DACQ) 
system. A brief description of the function of these subsystems 


is given along with block diagrams for each. 


A. TUBE AND PRESSURE VESSEL 


it. Resonant Tube 


The tube is a 1.22 m long, 3.81 cm inside diameter 
copper tube. The length was chosen to obtain one half of an 
acoustic wavelength - therefore a full thermal wavelength - at 
the fundamental. The tube is soft soldered to brass flanges 
which are attached to the pressure vessel at one end and a 
blanking plate at the other. The purpose of the blanking plate 
is to provide a rigid termination. An Endevco model 8510B-5 
piezoresistive transducer is threaded through the blanking plate 
and used to measure the acoustic pressure. A 1/8" diameter 
Stainless steel probe tube, which carries the TAC, also goes 
through the blanking plate. It is supported in the center of the 


plate by a rubber O-ring compression fitting. In order to allow 
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for elevated mean pressures within the tube, the volume behind 
the Endevco piezoresistive diaphragm is connected to the resonant 
tube via a capillary leak, which allows pressure equalization, 
but blocks acoustic pressure oscillations. The tube is wrapped 
with insulation material to ensure thermal isolation. 


(Figure 5). 


2. Pressure Vessel 


The ability to increase mean pressure required the 
acoustic driver (discussed below) to be enclosed in a pressure 
vessel. The vessel is made out of an aluminum cylinder joined to 
two aluminum end plates by threaded rods. The cylinder is 30.5 
cm long with a 25.4 cm inside diameter and 2.54 cm wall 
thickness. The end plates are 2.54 cm thick disks, 35.6 cm in 
diameter. The mean pressure inside the pressure vessel is sensed 
both by a dial pressure gauge and by an OMEGA Model PX304-150AV 
pressure transducer mounted in the back end plate. A 150 psi 


relief valve is mounted on the pressure vessel end plate. 


a. Driver 


The acoustic driver is a JBL Model 2445J compression 
driver with input provided by a Hewlett Packard Model 3314A 
Function Generator amplified by a Techron Model 7520 Power Supply 


Amplifier. The driver is mounted on the inside of the pressure 
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vessel end plate attached to the resonant tube. Provisions are 
made to enable mounting a positive displacement piston driver for 
work at high ambient pressures, though it was not used in this 


experiment. 


4. Plumbing 


Connections were made to pressurize the tube and 
pressure vessel, with provisions to equalize the pressure on both 
Sides of the acoustic driver. Connections were also made to 
attach a vacuum pump, with a needle valve for controlled 


backfilling, and also a system vent for outlet. 


B. TAC 


The TAC is a stack of 6 G-10 fiberglass plates, 1 cm long by 
2.54 cm wide (Figure 6). The center two plates are glued 
together with four type E (chromel-constantan) thermocouples in 
between. The thermocouples are connected in series to increase 
the output voltage and they allow the measurement of temperature 
difference across the TAC. The spacing between the center 
"Sandwich" and the remaining guard plates is 1.54 mm. A 
thermocouple is glued to the back edge of one of the guard 
plates, which allowed the measurement of the absolute TAC 
temperature. The center plate sandwich is 0.595 mm thick and the 


guard plates are 0.36 mm thick. 
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Thermoacous tic Couple (TAC) 


C. DACQ 


Data acquisition is centered around the HP 3457A multimeter 
(Figure 5). To reduce the number of connection lines in the 
figure, circled numbers represent connections to multimeter input 
channels and the circled "f" the front panel input. The Kikusui 


Model 6100 Oscilloscope was used to monitor system performance. 


The system reference temperature was read from a thermistor 
which was placed inside an aluminum block to ensure thermal 
stability (Figure 7). Junctions for tube and TAC temperature 
thermocouple wires are attached to this block, insuring the 


chromel-copper and constantan-copper junctions are isothermal. 


D. LINEAR POSITIONING 


Motive force for TAC positioning is provided by a Compumotor 
Model M83-135 Stepper Motor (25,000 steps/rev) (Figure 8) which 
is coupled to a Warner Electric 0.125 Lead Screw (8 turns/inch) 
on which rides a precision ball nut having zero lash. The ball 
nut is attached to an arm which slides on a support rod for 
Stability. The other end of the arm is attached to the stainless 
steel probe tube which holds the TAC. A Compumotor Model 2100 
Indexer serves as the interface between a Standard 286 (AT clone) 


computer and the motor. 


lig) 


09 


({BUIIIX9 


dWal Ff 


fies heres (£) dWaL OVL 


yooig wnuiwnypy 


oqnj sqoid 
ysnomy) uny 






LV 
rh cS 
077 a oD he, 

Ig 
preng 
LV OVE 


dWAL OVI 


Figure 7 


Temperature Reference 


18 


didD i. aNOWD LY 
UdsyAlTYd 











Jaxapu] 





OOTT [SPOW a 


AAAYOS AVAT 
Japoouq peondo LOAN 


T1IVd 


SEI-C8W 
Jojounduo-) 


9PNS 
La ae ee ee 


7X 
WUT] jvo1uryoopy AD 


\ 


poy uoddns 








A\DO 


Figure 8 


Linear Positioning System 


ee 


End of travel limit for both CW and CCW rotation and 


mechanical HOME (or reference) limit switches were tied directly 


to the indexer. All data runs start from the HOME position. 


TO 


more precisely define HOME, a 1000:1 optical encoder is used in 


conjunction with the mechanical home limit switch. Once the 
mechanical switch activates, rotation is continued until the 
appropriate index on the encoder wheel is found. Positional 
repeatability was consistently within 1 motor step, or 

approximately 0.127 microns of longitudinal travel, assuming 
mechanical switch was consistent enough to bring the optical 


encoder into the proper range. 


the 


Interface between the personal computer (PC), the indexer, 


and the other HP equipment is accomplished through a National 


Instruments IEEE 488 General Purpose Interface Bus (GPIB) card in 


the PC. 
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IV. MEASUREMENT 


Development and testing of the experimental apparatus with 
its complex control and data collection program was a major time 
consuming undertaking. However the result is that once the 
system iS manually pressurized, all subsequent operations are 


fully automated. 


Data acquisition is highly flexible. For example, the 
ability to adjust linear positioning step size ensured a data 
density sufficient to show the most detailed features of the 
temperature profile. Experimental results are presented in a 
serial manner with one feature changing per series. A summary 


figure leads to the discussion section. 


A. PROCEDURE 


Before data acquisition is started, the system is pumped down 
and backfilled with either helium or argon several times. 
Approximate resonance frequencies are determined next by driving 
the JBL with the noise source of the Signal Analyzer. The 
Endevco output is then fed back into the Signal Analyzer. The 
displayed signal peaks denote resonance frequencies. The driver 
input device is now changed to the function generator which is 


set to individual resonance frequencies to check for interference 
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between those resonances and higher harmonics of the drive signal 
(Figure 9). In general the harmonics were 18 or more dB below 


the fundamental. 


The gain on the voltage amplifier was set to a mid-range 
position and then left in that position for the duration of the 
data run. Determination of the magnitude of the amplification 


was unnecessary as acoustic pressure was measured directly. 


Automated data acquisition now begins. The program initially 
takes voltages from the ambient pressure transducer and directly 
converts them to pressure. The reference temperature is then 
read directly from the multimeter and converted to a effective 
reference voltage through a polynominal. The low order 
polynomial was numerically determined from NBS thermocouple 
tables over the limited temperature range of interest. Voltages 
from both the tube and TAC temperature thermocouples are read, 
added to the effective reference voltage, and converted to 
temperatures by an inverse polynominal. The drive voltage 
amplitude is set on the function generator and a search is made 
for resonance frequency by stepping on both sides of the 
previously determined natural frequency and reading the resultant 
acoustic pressure amplitude. The frequency producing the highest 
acoustic pressure amplitude is selected and set in the function 


generator. 
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Figure 9 


Drive Signal Spectrum 
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Once the preliminary readings are finished and the function 
generator set, the data run commences by reading the TAC 
thermocouple voltage and converting it to a temperature as 
discussed above. Then the TAC thermopile voltage is read and 


converted to a temperature differences as follows: 





Av dT (type E) 
At = NY 
4 dv Pacer 
Ar = TAC delta temp 
Av = TAC thermopile V 


dv TAC thermocouple voltage 

dT = TAC temperature 
The TAC is then moved to the next position, a wait time, obtained 
from Muzzerall [Ref. 3], activated and then readings are taken. 
Multiple data runs were made at constant static pressure and 
increasing acoustic pressure by resetting the indexer to the home 


position and incrementing the drive voltage on the function 


generator. 


B. RESULTS 


Experimental results are presented here in several series 
interspersed with theoretical predictions. The ordinate on all 
of the figures is the temperature adirterencecs| halos aan 
degrees Kelvin across the TAC. The abscissa is the product of 
the wave number and the distance of the center of the TAC from 
the rigid end of the resonance tube. All pressures, acoustic and 
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ambient are in kilopascals. The first plot (Figure 10) is the 
only one to show individual data points - all other data is shown 
by connecting lines. The purpose of this figure is to show the 
density of data points, which is similar for all the subsequent 


figures. 


The first series (Figures 11, 12, and 13) demonstrate the 
typical degradation of experimental results against theoretical 
prediction at increasing acoustic pressures. The next series 
(Figures 14 through 17) show the effect of shifting from the 
first through the third overtone. Remaining at the third 
overtone, the next series (Figures 17 through 21) show increasing 
ambient pressure in one atmosphere steps. The last series 
(Figures 22 through 26 ) shifts to Helium at one atmosphere and 


shows again increasing frequencies. 


C. DISCUSSION 


The graphical presentation of the "raw" experimental results 
show no sharp demarcation between agreement of theoretical 
prediction and measured data, but demonstrate a gradual degration 
of agreement at increasing frequencies and pressures. Several 


interesting features can be seen in the series. 


The first series demonstrate results akin those of Muzzerall 


(Ref. 3]. Specifically shown is good agreement at velocity nodes 
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Individual Data Points Plotted 
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Series la Increasing Acoustic Pressure 
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series lb Increasing Acoustic Pressure 
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Series lc Increasing Acoustic Pressure 
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Figure 14 


Series 2a Increasing Frequency 
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Figure 15 


Series 2b Increasing Frequency 
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Figure 16 


Series 2c Increasing Frequency 
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Figure 17 
Series 2d Increasing Frequency 


(AH) SBduauajysytTgG aunqeuaduas 


et °c 
aL4°2 
se°e™ 
66°T 
Te°t 
te°t 


#q4 veo’ 


Bd @8c 7H Ser? 


XA 





IWINSWITaSdxXS NOSaY 


(HW) SB0uUavayJTGQ aunjeusdual 


Figure 18 
Series 3a Increasing Static Pressure 
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Series 3b Increasing Static Pressure 
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Series 4b Increasing Frequency 
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(KX in multiples of n7l) and degradation at the pressure nodes 
(KX at (2n+1)(7T/2)). Similar skewing to what he showed was also 


evident. 


A second feature present in all of the series was an 


asymmetry in the magnitude of the maximum temperature 


differences. The magnitude of the positive temperature 
G@itference (T, - TT.) was always less than that of the negative 
eee 2, ) difference. 


The shape of the curves also exhibited interesting 
properties. The curve shape was remarkably similar from cycle to 
cycle on the same side of the Ar axis, but varied greatly from 


above to below the axis. 


Examination of the maximum temperature differences vs 
pressure amplitude (Figure 27) reveals the rate at which the 
temperature difference falls off as the acoustic and mean 
pressure increases. Several of the data sets start to show a 
Saturation of the effect at high acoustic and mean or static 


pressures. 
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Maximum Temperature Difference vs Pressure Amplitude 
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V. SUMMARY/RECOMMENDATIONS 


The main accomplishment of this thesis was to put together a 
system enabling computerized measurement of temperature gradients 
developed across short plates by thermoacoustic heat transport. 
Measurements were made of the gradients under different static 
and acoustic pressures, different frequencies, and different 
gases. Comparison with theoretical values showed results Similar 
to those found by Muzzerall and further showed a gradual 
degradation of agreement from low to high acoustic pressures. 

The thermoacoustic effect becomes saturated at the higher values 
of acoustic pressure amplitude. Also, theoretical values and 
measured data diverge most severely in the area of pressure 
nodes. The causes of both of these effects bear further 


investigation. 
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